Evolution of the single-mode squeezed vacuum state in amplitude dissipative channel 
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Using the way of deriving infinitive sum representation of density operator as a solution to the 
master equation describing the amphtude dissipative channel by virtue of the entangled state rep- 
resentation, we show manifestly how the initial density operator of a single-mode squeezed vac- 
uum state evolves into a definite mixed state which turns out to be a squeezed chaotic state with 
decreasing-squeezing. We investigate average photon number, photon statistics distributions for this 
state. 



I. INTRODUCTION 

Squeezed states are such for which the noise in one 
of the chosen pair of observables is reduced below the 
vacuum or ground-state noise level, at the expense of 
increased noise in the other observable. The squeezing 
effect indeed improves interferometric and spectroscopic 
measurements, so in the context of interferometric detec- 
tion and of gravitational waves the squeezed state is very 
useful [H, Q. In a very recently published paper, Agarwal 
[3^ revealed that a vortex state of a two-mode system can 
be generated from a squeezed vacuum by subtracting a 
photon, such a subtracting mechanism may happen in 
a quantum channel with amplitude damping. Usually, 
in nature every system is not isolated, dissipation or de- 
phasing usually happens when a system is immersed in 
a thermal environment, or a signal (a quantum state) 
passes through a quantum channel which is described by 
a master equation [4]. For example, when a pure state 
propagates in a medium, it inevitably interacts with it 
and evolves into a mixed state [5]. Dissipation or de- 
phasing will deteriorate the degree of nonclassicality of 
photon fields, so physicists pay much attention to it p-|8|. 

In this present work we investigate how an initial 
single-mode squeezed vacuum state evolves in an ampli- 
tude dissipative channel (ADC). When a system is de- 
scribed by its interaction with a channel with a large 
number of degrees of freedom, master equations are set 
up for a better understanding how quantum decoherence 
is processed to affect unitary character in the dissipation 
or gain of the system. In most cases people are inter- 
ested in the evolution of the variables associated with 
the system only. This requires us to obtain the equa- 
tions of motion for the system of interest only after trac- 
ing over the reservoir variables. A quantitative measure 
of nonclassicality of quantum fields is necessary for fur- 
ther investigating the system's dynamical behavior. For 
this channel, the associated loss mechanism in physical 
rocesses is governed by the following master equation 



dp {t) 
dt 



n (2apa^ — a^ap — pa^a) , (1) 



the rate of decay. We have solved this problem with use 
of the thermo entangled state representation [9]. Our 
questions are: What kind of mixed state does the initial 
squeezed state turns into? How does the photon statistics 
distributions varies in the ADC? 

Thus solving master equations is one of the fun- 
damental tasks in quantum optics. Usually peo- 
ple use various quasi-probability representations, such 
as P-representation, Q-representation, complex P- 
representation, and Wigner functions, etc. for converting 
the master equations of density operators into their corre- 
sponding c- number equations. Recently, a new approach 
p, lOl, using the thermal entangled state representation 
[l2| to convert operator master equations to their c- 
number equations is presented which can directly lead to 
the corresponding Kraus operators (the infinitive repre- 
sentation of evolved density operators) in many cases. 

The work is arranged as follows. In Sec. 2 by virtue 
of the entangled state representation we briefly review 
our way of deriving the infinitive sum representation of 
density operator as a solution of the master equation. 
In Sec. 3 we show that a pure squeezed vacuum state 
(with squeezing parameter A) will evolves into a mixed 
state (output state), whose exact form is derived, which 
turns out to be a squeezed chaotic state. We investigate 
average photon number, photon statistics distributions 
for this state. The probability of finding n photons in this 
mixed state is obtained which turns out to be a Legendre 
polynomial function relating to the squeezing parameter 
A and the decaying rate hz. In Sec. 4 we discuss the 
photon statistics distributions of the output state. In 
Sec. 5 and 6 we respectively discuss the Wigner function 
and tomogram of the output state. 



II. BRIEF REVIEW OF DEDUCING THE 
INFINITIVE SUM REPRESENTATION OF p (t) 

For solving the above master equation, in a recent 
review paper [l^ we have introduced a convenient ap- 
proach in which the two-mode entangled state [ijj, 



where p is the density operator of the system, and Hi is 
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is employed, where is a fictitious mode independent 
of the real mode , [a^a^] =0. \r] = 0) possesses the 
properties 



a\r] = 0) = a'^lf] = 0), 
a^r] = 0)=a\r] = 0), 
(ata)^|7? = 0) = (ata)^|7? = 0). 



(3) 



Acting the both sides of Eq.([T]) on the state \r] = 0) = \I) , 
and denoting \p) = p\I)^ we have 



d_ 
~dt 



\p) = K, (2apa^ — o)ap — po) a) \I) 

= n (2ad — a^a — a^a) \p) , (4) 

so its formal solution is 

\p) = exp [f^t (2aa — a^a — a^a)] |po) , (5) 

where |po) = Po |^) , Po is the initial density operator. 

Noticing that the operators in Eq.(j5j) obey the follow- 
ing commutative relation, 



[aa,a'''a] = [aa,a'''a] = aa 



and 



-, aa 



as well as using the operator identity 
(which is valid for [A, B] = rB)^ we have 



-2Kt 



where T = 1- e'^^'K 

Then substituting Eq.® into Eq.© yields [Tsl 



(6) 



(7) 



(8) 



(9) 



^ rnfn 



n ^ — Ktcv a 



n=0 



^ ^e-'^'^^a^poa^'^e-"'^^ |/) , (10) 



n=0 



which leads to the infinitive operator-sum representation 
of p. 



n=0 



where 



(11) 



(12) 



We can prove 
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2nKt , n . —2Kta^o 



= : e(^'^*-i)^^« : g-^'^^'^^^ = 1, (13) 

where : : stands for the normal ordering. Thus Mn is a 
kind of Kraus operator, and p in Eq. (pT]) is qualified to 
be a density operator, i.e.. 



Tr[p{t)] = Tr 



n=0 



Trpo. (14) 



Therefore, for any given initial state po, the density op- 
erator p{t) can be directly calculated from Eq. (pT]) . The 
entangled state representation provides us with an ele- 
gant way of deriving the infinitive sum representation of 
density operator as a solution of the master equation. 



III. EVOLVING OF AN INITIAL SINGLE-MODE 
SQUEEZED VACUUM STATE IN ADC 

It is seen from Eq. (pT]) that for any given initial state 
po, the density operator p{t) can be directly calculated. 
When po is a single- mode squeezed vacuum state. 



/ tanhA 



Po = sechAexp \^ — - — a^'^ j |0) (0| exp \^ ^ 



( tanh A 



(15) 



we see 



n=0 



X (0| exp 



/ tanh A 



V 2 



^2 \ ^\n^-nta^a^ 



\ 2 

Using the Baker-Hausdorff lemma [l5| , 

A^ 



e^^Be-^^ = B^X 

we have 

/ tanhA 



A,B 



2! 



A 5 



(16) 



. (17) 



a exp ^ — - — a 



t2 



|0) 



tanh A t2 tanh A t2 tanh A t2 



" |0) 



tanh A ^12 / x ^ 

= (a + aUanhA) |0).(18) 

Further employing the operator identity 



(/ia + i 



— a ^ i\ — a' 

V 2m y 

(19) 
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where H^{x) is the Hermite polynomial, we know 
(a + tanh A)^ 

• Hn i\l — — —a^i 



tanh A 



From Eq.([T8]), it follows that 



a^e^^ |0) 



2 tanh A 



tanh A 



tanh A 



a"^ bo) 



xHJ^,|^aA\0). (21) 



By comparing Eq. (p!5|) with (|23|) one can see that after 
going through the channel the initial squeezing parame- 
ter tanh A in Eq.([T5|) becomes to B= i^_y?2 tanh A ^ 

|0)(0| ^ 1 gataln(BT'tanhA) ^ ^^^^^-^ ^^^^^ 

— 2Kt 

(mixed state), due to > 0, we can prove tanh^ a ^ 
1, which means a squeezing-decreasing process. When 
f^t = 0, then = and B = tanh A, Eq.l^ becomes the 
initial squeezed vacuum state as expected. 

It is important to check: if Ti p{t) = 1. Using Eq. (|22]) 

and the completeness of coherent state J ^ (^| = 1 
as well as the following formula [18,] 



On the other hand, noting e '^^^ ^a^e" 



7T 



(26) 



ae « and the normally ordered ^^ose convergent condition is Re(C ± / ± 5) < 



form of the vacuum projector |0) (0| = : e ° - , we have (9^) < 0, we really 



see 



p (t) = sechA y — -e-'^^^^a^e 

n=0 



Trp(t) = W [ i£(^|ef°t%»^»i"Kt-hA)gfa^|^^ 



tanh A ^2 



a ]n —Kta^ a 



w 



— sechA 



(T'tanhA)'" e-2^^at2 ^anh a 



n=0 



= 1. 



(BT'tanhA- 1)^ -B^ 



(27) 



X Hn\-i 



— sechA 



tanh A 



|0) (0| 



tanh A 



-ae 



n=0 



(T^ tanh A) 
2^n! 



n p-2Kt 



(a2+at2)tanhA , 
-i^ 1 n'r. 



■ : e 



so /) (t) is qualified to be a mixed state, thus we see 
an initial pure squeezed vacuum state evolves into a 
squeezed chaotic state with decreasing-squeezing after 
passing through an amplitude dissipative channel. 



IV. AVERAGE PHOTON NUMBER 



X Hr> 



tanh A 



tanh A 



then using the following identity [ij 



2^n! 



(x^ + y^) - 2txy 



and e^^^^ =: e(^^ i)a^a. ^ finally obtain the expres- 
sion of the output state 

^ |ygfat2^ataln(BT'tanhA)gfa2^ ^24) 

with T' = 1- and 



sechA 



-2Kt 



tanh A 



Vl-T'^anh^A 1 - tanh^ A 



-Kt\ . Using the completeness relation of coherent state 
j and the normally ordering form of p {t) in Eq. 
(22) ([22]) . and using e^^ a^e~i^ = a''"+Ba, as well 

as e^'^^"K*^"h^)ate-^'^^"(^^'*^"h^)=atBr tanhA, we 
have 



Tr [p{t) a^a) 

J fl,e^{^"+^n (^1 (atT'tanhA + ^) e-'"'"^*-^^) \z) 
j ^ (|^|2r'tanhA + z2) 

(28) 



, (23) 



(25) 



= WQ 

= we 

X exp 



(6T'tanhA-l)|2|2 + ^(z*2 + 22^ 
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FIG. 1: (Color online) The average n{hit) as the function of 
for different values of squeezing parameter A (from bottom 
to top A = 0,0.1,0.3,0.5, 1.) 



In order to perform the integration, we reform Eq. ([28|) 
as 

Tr{p{t)a^a) = l^B tanh A^ + |^ } 

X y^^exp [(BrtanhA-l + /)|zp 



+ ^(.*^ + (l + .).^) 

1 - BT' tanh A 
(BT' tanh A - 1)^ -B^ 



(29) 



f=s=0 
- 1 



in the last step, we have used Eq. (j27j) . Using Eq.(j29]), we 
present the time evolution of the average photon num- 
ber in Fig. 1, from which we find that the average pho- 
ton number of the single-mode squeezed vacuum state in 
the amplitude damping channel reduces gradually to zero 
when decay time goes. 



V. PHOTON STATISTICS DISTRIBUTION 

Next, we shall derive the photon statistics distribu- 
tions of p{t). The photon number is given by p{n,t) = 
{n\p{t)\n). Noticing a"*"^ I n) = y/ {m -\- n)\ / n\ \m -\- n) 
and using the un- normalized coherent state \a) = 
exp[aat] |0), ^ ^ leading to |n) = \a) \a=o , 



{{f3 \a) 



, as well as the normal ordering form of 



p{t) in Eq. (j22j), the probability of finding n photons in 
the field is given by 



{n\p{t)\n) 
W 

W 

n! da^ 



I 62 



exp 



BT' tanh A + ^/3*2 ^ 



-a 
2 



2 



Note that 



f a"''aln(6T'tanhA) f a^l t ^ |at2 ln(BT' tanh A) f 



{n\ pit)\ny = {n\ p{t)Un) = {n\ pit)\n) 



Qn-\-n 

[n/2] 

2"n! 



I 



n! 



22^ {m\y{n-2m)\ 
we derive the compact form for p (n,t), i.e.. 



(31) 



p(n,t) 
W ( ^ 



2 ] da^ 

[n/2] 



(BrtanhA)"" ^ 



-2T' tanh A/3*q;-/3*^-q;^ 



n!(rtanhA)" 



^q22^ (m!)^(n-2m)! 



(32) 



Using the newly expression of Legendre polynomials 
found in Ref. ml 



[n/2] 



'g22™ (m!) (n-2m)! 



= Pn (x) , 



(33) 

we can formally recast Eq. (|32]) into the following compact 
form, i.e., 

p{n,t) = W (^e-'^V-Btanh Xy Pn (^e'^^TV-B tanh a) 



note that since V— Btanh A is pure imaginary, while 
p (n, t) is real, so we must still use the power-series expan- 
sion on the right-hand side of Eq. (|32]) to depict figures of 
the variation of p (n, t). In particular, when t = 0, Eq. (|32l 
) reduces to 

/ \ 1 A / 1 A\n T n!(i tanh A) 
p(n,0) = sechA(tanhA) lim > ^ 



22fcfc!fc! sechA tanh^^ A, n = 2k 







n = 2A: + l ' 



(34) 



which just correspond to the number distributions of the 
squeezed vacuum state [H, |23[. From Eq. (|34|) it is not 
difficult to see that the photocount distribution decreases 
as the squeezing parameter A increases. While for f<it 
oo, we see that p (n, oo) = 0. This indicates that there is 
no photon when a system interacting with a amplitude 
dissipative channel for enough long time, as expected. 
Fig. 2, the photon number distribution is shown for 
ifferent Hit. 



FIG. 2: (Color online) Photon number distribution of the 
squeezed vacuum state in amplitude damping channel for A = 
1, and different K,t: (a) K,t = 0, (b) nt = 0.5, (c) /^t = 1 and 
{d) Kt 2. 



FIG. 3: (Color online) Wigner function of the squeezed vac- 
uum state in amplitude damping channel for A = 1.0, different 
Kt: (a) Kt = 0.0, (b) nt = 0.5, (c) = 1, and {d) nt = 2. 



VI. WIGNER FUNCTIONS 

In this section, we shall use the normally ordering for 
of density operators to calculate the analytical expression 
of Wigner function. For a single-mode system, the WF 
is given by [M] 

W (a, a\t)= e^l-l' j ^ p {t) \p) e-^C/^"*-/^*"), 

(35) 

where \j3) is the coherent state [l^,|20| . From Eq.(|22) it 
is easy to see that once the normal ordered form of p {t) is 
known, we can conveniently obtain the Wigner function 
ofp{t). 

On substituting Eq.^ into Eq. (|35l) we obtain the WF 
of the single-mode squeezed state in the ADC, 



W{a,a*,t) 



-2 (/3a' 



(5* a) 
W 



■ exp 



(1 +6T'tanhA) 



3*2 



: exp 



. 2\a\ 

7r^/(l + BT' tanh A)^ -B^ 

-2 (1 + BT' tanh A) \a\^ + B (a*^ 



X exp 



(1 + BT' tanh A)^ - B^ 



-36) 



In particular, when t = and t oo, Eq. (|36|) 

reduces to VK(a,a*,0) = exp[-2 cosh2A + 

sinh2A], and W(a,a*,oo) = :^exp -2|a|^ , 

which are just the WF of the single- mode squeezed vac- 
uum state and the vacuum state, respectively. In Fig. 3, 
the WF of the single-mode squeezed vacuum state in the 
amplitude damping channel is shown for different decay 
time hvt. 



VII. TOMOGRAM 

As we know, once the probability distributions Pq (xq) 
of the quadrature amplitude are obtained, one can use the 
inverse Radon transformation familiar in tomographic 
imaging to obtain the WF and density matrix [25]. Thus 
the Radon transform of the WF is corresponding to the 
probability distributions Pq (xq). In this section we de- 
rive the tomogram of p {t). 

For a single-mode system, the Radon transform of WF, 
denoted as 1Z is defined by ^ 

^ = jHQ-fq'- 9P') Tr [A p {t)] dq'dp' 

= f,9 f,9 (l\ P (*)] =/.5 P l9) 

where the operator \q) ^ g ^ g{q\ is just the Radon trans- 
form of single- mode Wigner operator A (/3), and 



f,9 



A exp 



B 



2b' 



,t2 



|0): 



(38) 



as 



well as 5 = / 

-1/4 „,.„r „2 /o /i'2 , 



- ig, A 

[iT{f^g^)]-"%M-<lV2{f^9^)Y Thus the 
tomogram of a quantum state p (t) is just the quantum 
average of p(t) in \q)f^g representation (a kind of 
intermediate coordinate- momentum representation) [23| • 
Substituting Eqs.(j24|) and (jSH]) into Eq.(j3ll), and using 
the completeness relation of coherent state, we see that 
the Radom transform of WF of p {t) is given by 



f,g 



exp 



E 

2q^Q 
E\B\ 



[E\B\ 

(T tanh X + Q-QT^ tanh^ A) I , (39) 



6 



where we have used the formula 
exp[- \a\^ /2 - |7|^ /2 + a*7], as weh as 



E 




( BrtanhA)^ 



In particular, when t = 0, (T = 0), then Eq. (|39|) reduces 



B_ ^2i4 



A^sechA 



•tanhA| 



X exp 



(7^ (^^2 ^ 5*2 + 2 tanh A^ tanh A 



l + e^^^l^rtanhA 



which is a tomogram of single-mode squeezed vacuum 
state; while for hzt oo,(T = 1), then lZ{q)j ^ = A^, 
which is a Gaussian distribution corresponding to the 
vacuum state. 



In summary, using the way of deriving infinitive sum 
representation of density operator by virtue of the entan- 
gled state representation describing, we conclude that in 
the amplitude dissipative channel the initial density op- 
erator of a single-mode squeezed vacuum state evolves 
into a squeezed chaotic state with decreasing-squeezing. 
We investigate average photon number, photon statistics 
distributions, Wigner functions and tomogram for the 
output state. 
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